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Abstract;  A second-order  linearized  theory  of  two-dimensional  cavity 
flov7S  is  used  to  study  the  nonsteady  response  of  hydrofoil 
lift  and  moment  to  sinusoidal  and  sharp-edged  streamwise 
gusts.  These  gust  patterns  are  assumed  to  be  frozen  and 
to  be  convected  at  the  steady  free-stream  velocity.  The 
analysis  is  restricted  to  the  case  of  zero  cavitation  number 
corresponding  to  an  infinitely  long  cavity  in  an  otherwise 
unbounded  flov;.  For  sinusoidal  gusts,  lift  and  moment 
response  functions  are  presented  for  the  entire  range  of 
reduced  frequency.  In  addition,  transient  iift  and  moment 
responses  are  tabulated  for  the  reactions  which  occur  after 
the  foil  encounters  a sharp-edged  gust.  These  calculations 
are  carried  out  for  those  terms  in  the  solution  which  result 
from  the  nonsteady  do\-mwash  on  the  wetted  surface  of  the 
foil  which  are  due  to  the  direct  action  of  the  gust  on  the 
inclined  wetted  surface.  They  provide  a direct  cavity  flow 
analog  of  the  Horloc.k  function  of  airfoil  theory.  The 
present  study  cf  a cavity  flow  lias  revealed  an  added  non- 
steady effect  resulting  from  an  interaction  between  the 
gust  and  the  cavity,  Tlicoretirnl  results  for  this  part  of 
the  solution  arc.  iiresented  in  this  report  which  show  a 
response  at  Lv/<ca  the  input  frequc--ncy  of  a sinusoidal  gust. 
Kumcricnl  evaluation  of  this  part  of  the  theoiy  has  not  been 
carried  out;. 
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INTRODUCTION 


The  present  study  Is  motivated  by  recent  demonstrations  that  accurate 
analyses  of  nonsteady  blade  forces  caused  by  inflow  distortions  upstream  of 
propellers  [l,2]*and  axial-flow  turbomachines  [3,4]  depend  upon  the  stream- 
wise,  as  well  as  the  transverse,  component  of  gust  velocity  relative  to  tlie 
blade.  In  these  experimental  and  theoretical  studies,  which  Involve  non- 
cavitating  flows,  the  inclusion  of  the  streamwise  gust  component  has  been 
found  essential  for  prediction  of  nonsteady  force  trends  with  advance  ratio. 
Another  reason  for  the  present  analysis  is  the  current  Interest  in  high-speed 
hydrofoil  boats  with  supercavitating  foils.  In  such  craft,  the  <"orces  on 
supercavitatlng  propeller  blades  or  ocean  wave  effects  on  the  foil  loads 
might  also  depend  on  the  stream\<7ise,  as  well  as  the  transverse,  component  of 
gust  velocity.  Finally,  the  possibility  that  streamwise  gusts  might  also  be 
generated  in  the  course  of  nonsteady  cavity-flow  experiments  in  water  tunnels 
is  another  reason  for  considering  this  problem. 

The  present  study  is  only  a start  in  the  analysis  of  the  nonsteady ' cavity 
flow  problems  noted  above.  With  respect  to  the  first,  the  effects  of 
neighboring  blades  are  not  considered,  with  respect  to  the  second,  the  presence 
of  the  free  surface  is  neglected,  and  with  respect  to  the  third,  the  effect  of 
tunnel  walls  has  not  been  Included,  The  aim  here  is  to  obtain  the  cavity-flow 
analog  of  the  Horlock  function  [5]  for  nonsteady  cavity  flows,  and  from  this 
new  result  to  derive  nonsteady  response  functions  for  lift  and  moment  to 
sharp-edged  gusts.  The  entire  analysis  is  restricted  to  two-dimensional  flows 
at  zero  cavitation  number.  In  this  case,  the  cavity  length  is  Infinite.  The 
linearized  theory  of  cavity  flows  is  used  to  obtain  the  desired  results. 

As  is  also  true  of  Horlock' s analysis  [5],  the  nonsteady  response  of 
streamwise  gusts  in  cavity  flows  is  a second-order  effect.  Therefore,  the 
present  analysis  centers  around  the  systematic  development  of  first,  and  then 
second,  order  approximations.  As  is  convenient  in  linearized  analyses  of  sucli 
nonsteady  cavity  flows,  the  acceleration  potential  is  employed.  Tlie  method  of 
the  present  study  follows  an  analogous  treatment  of  nonsteady  airfoil  theory 
[6].  However,  the  presence  of  the  long  cavity  in  the  hydrofoil  case  brings 
certain  aspects  to  the  theory  which  require  that  we  reformulate  the  analysis 
from  the  beginning. 

FLOW  GEOMETRY  AND  BAS  .C  EQUATIONS 

We  shall  consider  a two-dimensional  flat-plate  hydrofoil  of  chord  length  s 
at  an  attack  angle  £0,  v;itli  respect  to  tlie  free-stream  flow  direction.  The  cavity 
springs  from  the  nose  and  tail  of  the  profile,  and  it  extends  to  infinity  in  the 
downstream  direction.  Only  the  underside  of  the  foil  is  wotted  by  the  flow,  and 
this  surface  lies  along  the  x axis.  Figure  1 shows  the  general  features  cif  the 
flow  geometry  v;ith  respect  to  x-y  coordinates.  The  free-stream  velocity,  , 

in  Figure  1.  includes  the  sinusoidal  gust,  Uw(x,t)  , which  is  supc^rimposed  on 
the  free-stream  velocity  U.  The  gust  velocity  is  the  only  time-dependent  input 
disturbance  to  an  otherwise  steady  flow.  The  attack  angle  f;a  is  fixed,  and  the 
pressure  inside  the  cavity  is  constant  aiul  equal  to  the  free-stream  static 
pressure,  P^.  The  quantity  C.  cliaractcti  zes  the  magnitude  of  all  disturbc.nces 
imposed  upon  the  free-stream  flow. 

The  velocity  at  any  point  in  the  flov:,  q(x,y,t)  , has  two  coi.iponcni s 

* 

Numbers  in  pareuiLeses  refer  to  documoiits  given  In  tlie  list  of  references. 
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q = U (cosea  + w + u) 

X X 

and  (1) 

= U(sinea  + + v)  . 

Each  of  these  has  a steady  free-streara  part,  the  periodic  gust,  and  a disturb- 
ance component  created  by  the  foil  and  cavity.  The  disturbances  u and  v 
vanish  far  from  the  foil.  Moreover,  the  static  pressure  at  any  point  in  the 
flow  is  P(x,y,t).  It  can  be  represented  as  the  sum  of  the  constant  free-stream 
pressure  Pg,  and  a disturbance  pressure  p(x,y,t)  created  by  the  foil.  Thus, 


P(x,y,t)  = Pq  + P(x,y,t) 


(2) 


The  disturbance  p(x,y,t)  also  vanishes  at  points  far  from  the  foil  and  on  the 
cavity. 

For  an  incompressible  fluid  having  density  p,  conservation  of  momentum 
requires  that 


9q  , VP 

Tt  + = " T 


(3) 


Conservation  of  mass  provides  that 


div  q = 0 


(4) 


If  the  flow  is  Irrotational, 


curl  q = 0 


(5) 


In  a two-dimensional  flovj.  Equations  (4)  and  (5)  imply  that  a gust  in  the 
x-dlrectlon  should  be  associated  with  one  in  the  y-direction.  In  nonsteady 
airfoil  theory,  the  response  to  a constant  amplitude  vertical  sinusoidal  gust 
is  given  by  the  well-known  Soars'  function  [7],  The  response  to  a similar 
streamwise  gust  has  been  determined  by  llorlock  [5],  with  a further  extension 
and  interpretation  given  by  Morfey  [8].  Commerford  and  Carta  [9]  have  studied 
more  general  two-dimensional  sinusoidal  transverse  and  chordwise  gusts  when 
the  conditions  of  irrotationality  and  continuity  arc  explicitly  satisfied. 

T!iey  find  that  gusts  of  constant  amplitude,  independent  of  y,  are  not  strictly 
pennisslble.  Both  components  should  show  tranr:versc  amplitude  variations 
represented  by  multiplicative  factors  of  the  form  exp(luy/l)).  They  also  find 
that  the  phase  of  the  liorizontal  gust  leads  that  of  the  vertical  gust  by  90 
degrees.  The  neglect  of  those  conditions  in  analyses  by  Scars,  Horlock  and 
others  appears  to  be  justified  by  the  fact  that  small  voriicity  in  tiie  flow 


- 5 - 


June  4,  1974 
BRPrlhm 


does  not  affect  the  lift  on  a thin  airfoil.  Evidently,  independent  analyses 
of  constant  amplitude  vertical  and  horizontal  gust  response  can  provide 
practically  meaningful  results.  Appropriate  combinations  of  such  transverse 
and  horizontal  gust  response  functions  can  then  be  used  to  approximate  the 
response  to  practical  situations  [5]. 

PERTURBATION  EQUATIONS 

The  effect  of  the  gust  will  be  studied  under  the  assumption  that  the 
hydrofoil  and  its  cavity  are  equivalent  to  a slender  body  so  that  disturbances 
created  by  the  foil  will  be  much  smaller  than  the  free-stream  velocity  U.  The 
parameter  e will  be  used  to  characterize  the  magnitude  of  all  disturbances. 

In  the  case  of  the  gust,  the  chordwise  component  which  is  parallel  to  the 
wetted  surface  of  the  flat  plate  foil  is  the  interesting  quantity.  We  will 
specify  it  by 


w 


= ewj^  = 


eu  e 
o 


30J(t  - 


(6) 


where  the  constant  Uq  Is  the  chordwise  gust  amplitude,  to  is  the  circular 
frequency,  t is  the  time  and  j = ■'^-1.  It  follows  that 


w = ew,  tanea  = e aw,  + . . . 
y 1 1 


(7) 


It  may  be  remarked  that  if  Wy.  is  not  at  least  of  first  order  in  £,  it  follows 
from  the  continuity  condition  that  it  must  be  zero,  which  is  contrary  to 
hypothesis. 

Now,  let  disturbances  created  by  the  foil  be  represented  as 


u = euj^  + e u^  + . . . , 


V = cv^  + c v^  + . . . , 


(8) 


and 


p = epj^  + c p^  + 


/ 


Expanding  the  trrgnometric  functions  in  Equations  (1)  in  pov.’ors  of  £ and  using 
Equations  (6),  (7)  and  (8)  in  Equations  (I),  one  gels  the  component  velocities. 


and 


q^  = U fl  + C(uj^  + wp  + (U2  - J-)  + ...]  , 


q = U [t  (a  + V,)  i c.  (v  + aw  ) + 

y J Z i 


(9) 
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Now,  let  us  substitute  Equations  (9)  into  the  irrotationality  condition  of 
Equation  (5).  Then,  group  all  terms  Recording  to  the  power  of  e multiplying 
the  terms  of  each  group.  Since  curl  q = 0 for  all  admissible  values  of  e, 
it  follows  that 


and 


0 


3x 


(v^  + oew^) 


0 


J 


(10) 


Similarly,  the  continuity  equation.  Equation  (4),  becomes 


and 


3 

_ (^1  + Wf)  + -^  = 0 


3v 

— - + — - = 0 
3x  3y 


(11) 


The  same  process  applied  to  the  equations  of  motion  also  leads  to  pairs  of 
first  order  and  second  order  equations. 


First  Order  Equations 


U 3t 


\ 


/ 


(12) 


Tlic  gust  velocity,  because  of  its  form,  has  zero  linearized  total  derivative, 
and  it  does  not  appear  in  Equations  (12).  Moreover,  v/h^n  tlie  first  of 
Equations  (11)  is  combined  with  the  sum  of  Equations  (12),  after  the  first  is 
differentiated  w’ th  respect  to  x,  and  the  second  v;ith  respect  to  y,  one  finds 
that 


(i  <L+  3_] 

U 3t  3x 


3x 


0 


/ 
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again,  because  = Uq  e ^ . Therefore,  we  can  define  a harmonic 

function  such  that  the  gradient  of  this  acceleration  potential  gives  the 
disturbance  acceleration  created  by  the  foil: 
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Again,  differentiation  of  the  first  member  of  this  pair  with  respect  to  x, 
and  the  second  with  respect  to  y and  use  of  the  second  order  continuity 
equation,  leads  to 

n P,  (w,  + u )^  + (v,  + a)^ 

V < 5 > = 0 

pU 


As  before,  we  set 


2 2 

?2  (”2. 

+ _ 

pU 


V 

The  integral  of  this  equation  is 


pU 


(w^  + Uj^)^  + (v,  + a)^ 


When  Uq  = 0 the  flow  must  become  steady  and  so  f(t)  is  a constant  which  we 
will  put  equal  to  zero.  At  upstream  infinity,  P2  = = v-j^  = 0,  and  since 

the  value  of  <^2  is  at  our  disposal  there,  we  will  put  “ “ (w^+a^)/2.  As 

will  be  seen  below,  this  choice  of  insures  that  P2  is  zero  on  the 

trailing  edge  of  the  foil,  and  of  greater  importance,  that  4>2  is  also  zero 
at  the  trailing  edge.  Had  we  chosen  C()2(“,t)  = 0,  the  value  of  po  would  still 
be  zero  at  the  trailing  edge,  but  (fiy  would  have,  been  equal  to  a^/2  tliere.  The 


present  choice  is,  therefore,  the  more  convenient  one. 
atlons  it  follows  that 


From  these  conslder- 


<^2  = - 


(Wi  + u^)^  + (vj^  + a)^ 


In  Equation  (17)  the  teim  w^  is  bounded  but  not  well  behaved  at  z = How- 

ever, it  loads  to  no  second-order  acceleration  components  anywhere  in  the 
flov7.  Its  effe?ct  will  be  treated  separately  below’. 

2 

Clcariy,  V (J>2  = 0 and  the  acceleration  components  can  be  expressed  in 
terms  of  the  Cauohy-Rieniann  equations  as 


9'J'2  ^’^''2 


X,,  8x  3y 


9<f>2  ^^P2 

9y  8x 
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As  was  the  case  for  the  first  order  system,  one  can  define  a complex 
acceleration  potential 

In  the  limit  of  steady  flow,  the  first  of  Equations  (16)  becomes 

8x  9x 


and  If  one  Integrates  from  z = -<»  to  z = x,  he  gets 

2 

steady  ^^2  2 


The  second  of  Equations  (16)  becomes 


'^2  steady  ^^2 


BOUKDAR/  CONDITIONS 

The  solution  of  the  present  problem  depends  upon  the  determination  of 
an  analytic  function 

F(z,t;c)  = eF^  + £^^2  + •••  » (IS) 

which  at  each  Instant  Is  regular  outside  the  foil  and  the  cavity.  The  complex 
potential  F will  satlsfj'  certain  bo<mdary  conditions,  which  In  keeping  with  thin- 
body  approximations  are  transferred  to  a cut  along  the  positive  real  axis  of  the 
z plane.  This  rdathematical  representation  of  the  foil  and  the  cavity  is  1] lus- 
tra ted  in  Figure  2. 

As  noted  previously,  the  wetted  surface  of  the  foil  is  positioned  on  the 
real  axis  (0  < x < s,  y - 0-).  V.'e  shall  denote  the  cavity  contour  by 
y = eC, (x),  where  and  C_  refer  to  the  upper  or  the  lower  contours,  respec- 
tively. The  boundary  conditions  which  apply  in  the  z plane  are: 

(j)  p(x.y.t)  = 0 for  0 < y = ef’+(x);  s < X,  y = eC_(x). 

(ii)  F(z)  continuous  at  x = s,  y = 0-  (Kutta  Condition). 


4 
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2 2 2 

(iii)  F(z)  = - e (w^  + a )/2  at  z = -°°  [because  of  Equation  (17)] 

(Iv)  qy  = 0 for  0<x<s,y=0-  . 

Finally,  if  Vf  is  the  downwash  on  the  wetted  surface  of  the  foil,  the 
vertical  acceleration  of  the  fluid  there,  is 

It  ^f  = ^ • 

These  boundary  conditions  can  now  be  separated  into  first  and  second 
order  systems  with  the  help  of  Equations  (8),  (9),  (11),  (17),  and  (13).  In 
order  to  carry  out  this  separation  for  (i) , the  conditions  must  be  transferred 
from  the  cavity  contours  to  the  top  or  bottom  of  the  cut  along  the  real  axis. 
This  step  is  carried  out  by  using  a MacLaurin's  series  in  y = £C^  : 

P^  = Pj^  (x,0+,t)  + £C^(x)  -|^  p^  (x,0-,t)  + ...  , n = 1,2  . 

The  first  and  second  order  boundary  conditions  take  the  following  forms. 


First  Order 


I(i)  Re  {F^(x,0,t)}  = 0,  0 < X,  y = 0+  ; E < X,  y = 0-  . 

l(ii)  Fj^(z,t)  continuous  at  x = s,  y = 0-  (Kutta  Condition) 

I (iii)  F^(-oo,t)  = 0 . 

I(iv)  Vj^(x,0-,t)  = -a  , 0 5 X < s 

Evidentally,  these  conditions  indicate  a steady  flow  for  the  first  ordei 
solution,  so  that  Fj^  = or  = uj^(x,y),  ~ ~V,(x,y).  It  is  to  be  noted 

that  1 (iv)  guarantees  the  satisfaction  of  condi tioir  (v)  above. 


Secotid  Order 


3v^  w^  + (Vj  + a)^ 

Fx  2 


II  (ii) 


V^(z,l)  continuous  at  x = s,  y = 0-  . 


i 

i 
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Il(lii) 

II  (iv) 

II  (v)  Im 


2 

— at  z = -«o 
< s,  y = 0-  . 
X < s,  y = 0- 
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! In  the  first  of  these  conditions,  the  first  order  result  = ~ V]^  and  cavity 

I pressure^condition  = 0 have  been  used.  As  a result,  II(i)  becomes 

(})2  = - w^/2  on  the  lower  surface  of  the  cavity  at  the  trailing  edge  because 
V]^  = - a and  C “ 0 there.  But,  since  ^2  is  continuous  at  tlie  trailing  edge, 
it  follows  from  Equation  (17)  that  P2  = 0 there.  The  foregoing  basic 
equations  and  boundary  conditions  permit  us  to  attack  the  first  and  second 
order  problems. 

I FIRST  ORDER  SOLUTION 

I 

j Although  the  first  order  solution  is  a well-kno\'m  result  of  linearized 

j cavity  flow  theory,  it  will  be  reviewed  in  this  section  because  the  second 

order  solution  requires  specific  first  order  results  and  the  theoretical 
approach  to  the  first  order  problem  is  also  useful  in  the  second  order  case. 

^ In  the  following,  a sequence  of  conformal  mappings  will  be  used.  The  region 

outside  the  wetted  surface  of  the  foil  and  the  cavity,  represented  by  the 
' cut  along  the  positive  real  axis  in  the  z-plane.  Figure  2,  is  transformed 

into  the  region  above  the  real  axis  and  the  unit  semicircle  about  the  origin 
in  the  ^-plane.  This  result  is  achieved  by  first  mapping  the  region  outside 
the  cut  in  the  z-plane  into  the  upper  half  of  the  V-plane  by  means  of  a 
square-root  relationsliip.  Then  that  part  of  the  real  V-axis  corresponding 
to  the  wetted  surface  of  the  foil  is  niappcd  into  the  upper  half  of  the  unit 
circle  in  the  ^-pi.ane  by  means  of  a Joukowski  transformation.  Thus,  one  can 
write  for  both  first  and  second  order  flol^7S, 


\/f-  V 4 (C  +i)  -i  . (19) 


The  V-  and  r;~planes  with  corresponding  points  on  the  cavity  and  foil  are  shouTi 
in  Figures  3 and  4,  In  these  mappings  the  complex  acceleration  potential 
F(z,t;c)  is  invariant  at  corresponding  points. 

An  already  noted,  the  first  order  flow  is  steady  and  “ tij  - ivj . 

In  the  r^-plane,  an  analytic  function  which  sati.sfics  the  boundary  conditions 
1(1)  througli  l(iv)  can  be.  found  by  inspection  This  function  is 


I 


It  can  be  seen  that  f is  purely  imaginary  v.’ticii  5 ‘t'ly  point  on  tlie  re.al  axis. 
Eviduiitally,  f satinfic:s  1(1).  Tills  function  has.  a leading,  edge,  singulari’y 


1 
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which  provides  for  the  branching  of  the  flow  at  the  nose  of  the  hydrofoil. 

At  C = 1>  changes  in  velocity  must  be  of  the  same  order  as  the  free-stream 
velocity  so  that  the  small-disturbance  assumptions  of  the  theory  break  down 
in  the  neighborhood  of  the  nose.  On  the  other  hand,  it  is  certainly  clear 
that  the  Kutta  condition  I(ii)  is  satisfied  at  the  trailing  edge,  ? = - 1. 

In  fact.  Re  f = 0 there.  Moreover,  as  can  be  seen  from  Equation  (19),  large 
values  of  |^|  correspond  to  large  values  of  |z|.  Therefore  I(iii)  is 
certainly  satisfied.  Finally,  at  points  on  the  wetted  surface  of  the  foil 
? = e^®  and  2f  = cot0/2-i.  Accordingly,  f can  also  be  made  to  satisfy  I(iv) 
if  we  put  f equal  to 


Fi  - U^(.)  . ^ (20) 


I 


One  can-use  Equation  (19)  for  corresponding  points  on  the  profile,  for 
which  ^ = e^  and  "J  z/s  = - ~\j  x/s,  to  show  that  cosO  = 1-  2"^  x/s.  It  then 
follows  from  Equation  (20)  that 


II 
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One  can  also  write  specific  expressions  for  the  cavity  contours  C^(x). 
However,  these  expressions  will  not  be  needed  for  the  nonsteady  solution, 
and  we  can  omit  this  detail. 

DISCUSSION  OF  SECOND  ORDER  BOUNDARY  CONDITIONS 

The  boundary  conditions  for  tlie  second  order  solution  as  defined  by  the 
five  conditions  II  can  be  split  into  uncoupled  steady  and  nonsteady  parts. 
For  the  steady  part  we  have: 


Steady  Boundary  Conditions 
^2  9v^  (v^  + a) 2 

(i)  ^2  2 " 9x  2 ’ ^ 


(ii)  W2(z)  continuous  at  x = s,  y = 0- 


(iii)  <}>2  = - — at  z = -00  (and  U2  = V2  = 0 there). 


(iv)  V2  = 0,  0 < X < s,  y = 0-  . 


If  allowance  is  made  for  the  somewhat  different  notation  used  here,  and  the 
fact  that  we  consider  the  special  case  of  the  flat  plate  profile,  it  can  be 
seen  that  these  conditions  correspond  to  those  given  by  Chen  [10].  In  view 
of  Chen's  results,  we  need  not  consider  the  steady  problem.  We  can  proceed 
directly  to  the  study  of  nonsteady  gusts. 


Nonsteady  Boundary  Conditions 


w;  0 

(i)  Re  {F2(z,t)>  “ " 5 ^ 


(ii)  F2(z,t)  continuous  at  x = s,  y = 0-  . 


(iii)  a = a = 0atx=-«> 
X y 


(iv)  V2  = -a  0 < X < s,  y = 0-  . 


(v)  im  j ” ^ ^ y ” 0-  . 

Innpertlon  of  the  five  nonsleady  boundary  conditions  listed  above  reveals 
that  (i)  and  (iv)  have  different  beliavinr  witli  respect  to  tlie  gust  frequeney 
as  indicated  by  t lie  p(>wcr  of  W|  viiicli  occurs  in  encii  one  of  tb.ese  expressions. 
In  particular,  (i)  shov;s  tiie.  douivied  fr.  qnency  dependence  of  2tu  ami  (iv) 
depends  on  O)  only.  In  ordei'  to  lianJlc  this  situation  ve  sliall  let  tlie  solution 
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be  the  sum  of  two  parts:  F2(z;u))  + F2(z;2w).  The  first  part  accounts  for 
the  nonsteady  downwash  generated  on  the  wetted  surface  of  the  foil  and  the 
second  part  concerns  the  Interaction  between  gust  and  cavity.  Then  we  can 
split  the  nonsteady  boundary  conditions  by  writing  the  boundary  conditions 
for  F2(z;o))  so  that  its  condition  (i)  reads 


Re  {F„  (z;o))  } = 0 ; x > 


' 

0 

• . y = ' 

0+ 

s 

0- 

. ✓ 

with  all  others  remaining  as  written  above.  For  the  case  of  F2(z;2oj)  we  will 
write  condition  (iv)  to  read 


I 

! 

i 


I 


1 


V2=0;0<x<s,y=0- 


with  all  others  remaining  as  written  above.  In  this  latter  case  we  note  that 
if  this  homogeneous  streamline  condition  is  satisfied  on  the  wetted  surface, 
then  the  homogeneous  condition  (v)  for  the  acceleration,  Im  {dF2 (z ; 2cJ) /dz}  = 0 
will  be  satisfied  identically  at  points  on  the  foil.  Hov.’ever  if  the 

acceleration  condition  is  satisfied  first,  one  must  still  see  to  it  that  the 

streamline  condition  on  V2  is  satisfied  explicitly. 

Vie  will  nov;  turn  to  the  dcteirmination  of  the  function  F2(z;w),  v;hich 

accounts  for  the  nonstcady  downv.'ash  on  the  foil  and  wiiich,  as  we  shall  see, 

loads  to  the  direct  cavity~flow  analog  of  the  Horlock  function,  T(oj). 

THE  ANALOG  OF  HOiy,OCK’S  FUNCTION 


As  we  have  noted  above,  the  boundary  conditions  for  FoCz;^))  are; 
(i)  Re  {F^Czfoj)}  = 0 ; X > 

(ii)  F.j(z;oj)  continuous  at  x = s,  y = 0- 


r 'i 

r 

0] 

0+ 

7 . y = < 

Ui 

0- 

1 j 

(iii)  = 0 


(iv)  V2  = -Cl  Wj^  ; 0 < X < s , y =■  0- 


(v)  Im  {dF2 (z ;w) /dz}  =0;0<x<s,y=0- 


llaviiig  "suppressed"  the  term  Wj  in  tliese  conditions,  \;c  niu5;t  also  suppress  the 
corresponding  tcrii  in  IqiiUioii  (17).  Thus,  for  this  part  of  the  nonsleady 
solution  Fijuatiun  (17)  is  repl, iced  hy 
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*2 


p uj  + (Vj  ^■a)\ 

77" 5 -Vl 


(17a) 


Thus,  ({>2  for  this  case  will  show  no  frequency  doubling. 

Comparison  of  these  nonsteady  second  order  conditions  on  F2(z;(d)  with 
the  first  order  conditions  considered  previously  suggests  that  the  function  f, 
introduced  for  the  first  order  solution,  should  be  adaptable  to  the  nonsteady 
problem.  Condition  (iv)  will  lead  to  some  differences  in  the  method  of 
solution,  but  it  is  clear  that  the  present  conditions  (i) , (ii),  and  (iii)  are 
satisfied  by  f.  It  can  also  be  seen  that  condition  (v)  is  also  satisfied  by 
f.  One  can  verify  this  by  using_the  fact  that  in  the  ^-plane,  the  real  and 
Imaginary  parts  of  the  product  e^®df/d?  give  the  radial  and  transverse 
components  of  acceleration.  By  using  this  relationship,  one  can  verify  that  f 
produces  no  radial  component  of  acceleration  on  the  unit  circle.  Therefore, 
this  function  produces  no  acceleration  component  on  the  wetted  surface  in 
the  E-plane,  except  possibly  at  the  leading  edge  where  f is  singular.  In  order 
to  adapt  f to  the  present  situation,  we  will  introduce  a function 


A(t) 


A e 
o 


jwt 


(24) 


The  factor  Aq (to)  will  be  determined  by  condition  (iv).  The  time  dependence  is 
separated  out  as  indicated  in  Equation  (24)  because  W] , as  defined  in  Equation 
(6),  can  also  be  separated  in  this  way.  Thus,  we  shall  write  the  nonsteady 
solution  as 


iA  ej<"' 


(25) 


In  order  to  determine  the  quantity  Aq  and,  thus,  complete  the  solution  of 
Equation  (25),  it  is  necessai'y  to  satisfy  nonstcady  condition  (iv)  by  integrating 
the  second  of  Eqiuations  (16)  . Again,  we  separate  thej  spatial  and  temporal 
factors  by  writing  V2  = v^(x,y)e-l^‘  and  11^2  ~ (xjiOc-i^*"  to  obtain 


Dv  9t|j 

.03  , o o 

j — V + — 

U o ox  3x 


(26) 


This  equation  can  now  be  integrated  from 
to  a point  near  or  on  the-  wctltai  .surface 
frequency. 


X -■  along  the  negative  real  axis 
of  the  foil.  Introducing  tiie  reduced 


* 

Note 


that  now  the 


term  w^  if;  supprci 


■sud 


indicated  in  r,(|uation  (17a). 
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, _ CJS 


based  on  the  entire  profile  chord  length,  one  can  write 


V = - e 
o 


-jk^  f dip  jk-^ 


(28) 


for  the  Integral  of  Equation  (26).  Assuming  that  this  integral  is  convergent, 
one  can  consider  x to  be  a point  on  the  wetted  surface  and  he  can  differentiate 
both  sides  of  Equation  (28)  with  respect  to  x.  But  on  the  foil,  nonstcady 
condition  (v)  indicates  that  dip^/d-^  = 0 there.  Therefore,  it  also  follows  that 
9vq/3x  = - (jk/s)vQ  on  the  wetted  surface.  Moreover  condition  (iv)  indicates 
that  Vq  = - aug  exp(-jkx/s)  when  0 < x < s,  y = 0-.  Therefore  Equation  (28)  is 
identically  satisfied  at  all  points  on  the  wetted  surface  of  the  profile.  This 
fact  permits  us  to  take  x = Of  in  Equation  (28)  without  loss  of  generality. 
However,  in  order  to  secure  convergence  in  the  evaluation  of  Equation  (28),  we 
integrate  by  parts.  We  take  the  principal  vai.ue  of  the  integral  by  indenting 
the  contour  of  integration  by  means  of  a small  semicircle  about  the  nose  in 
the  lower  half  z-plane.  Thus,  the  path  of  integration  goes  from  to  ee^^ 
and  then  around  the  semicircle  from  ce^^  to  Ee^2TT.  Finally,  one  takes  the 
limit  £ ->  0.  Thus,  Equation  (28)  becomes 


+ ik  f 
0- 


•1  ^ 
jk~ 

- au  - ip  (0+)  + jk  / e ^ ij;  (— ) d (— ) 
o ’^o'  / o s s 


(29) 


From  Equation  (25)  one  finds  that  ip^(0+)  ~ - Aq/2  and  tliat 


Now,  one  can  put  x/s  = - p and  substitute  the  above  relationships  into 
Equation  (29)  to  obtain 


- an  - — ° '^^1  + j k 
o 2 I 


-V  /~V~1  +'l/U  -f  1 _ . 

V 2 


“ikn  , 
c dn  7 


Q,  one  finds  that 


- 2cx  u 

o 

1 I jk  Ij  (k) 


Solving  for 
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where  the  integral  Is  of  obvious  definition.  The  close  resemblence  of 
Equation  (30)  to  Martin's  [11]  result  for  the  transverse  sinusoidal  gust 
is  noteworthy.  In  Reference  [11],  Martin  has  put 


L(k)  = 


1 + jk  I, 


and  noted  that  L(k)  is  the  cavity-flow  analog  of  the  Sears'  function.  Numerical 
values  of  L(k)  are  tabulated  in  Reference  12.  For  values  of  k < 20,  the 
integral  was  evaluated  by  numerical  integration  [13].  For  calculations 
involving  reduced  frequencies  exceeding  20,  the  quantity  1 + jklj^  was  approx- 
imated by  an  asymptotic  expansion  [13]. 

NONSTEADY  FORCE  AND  MOMENT  FOR  THE  FIRST  PART  OF  THE  SOLUTION 

The  nonsteady  hydrodynamic  force  and  moment  depend  on  the  quantity  <j)2 
for  points  on  the  wetted  surface  of  the  hydrofoil.  We  can  separate  the  time 
dependence  by  writing  (j'2  = have  done  with  other  second  order 

quantities.  One  can  use  Equations  (25),  (30)  and  (31)  to  express  (})q  on  the 
wetted  surface  as 


[)  = -au  - 1 . 

'^o  o'  V V X 


From  Equation  (17a)  one  can  write  the  pressure  coefficient  on  the  foil  as 


% = ^^>2  ■'-  2- Vl 


This  exinresslon  can  be  decomposed  into  a steady  part  Cp„  and  a nonsteady  part 
Cp^  The  steady  part. 


r - 9 

S'2  " ^2  " 2~  2“ 


would  be  used  for  the  steady  second-order  problem.  The  nonsteady  part, 


C =-2  (f*  +u,uc 

p o J o 


is  applicable  to  the  present  nonstendy  probleri.  In  particular, 
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C = + 2au 

Po 


L(k)  + e 


The  unsteady  section  lift  and  moment  coefficients  are  obtained  by  integrating 
this  pressure  over  the  wetted  surface  of  the  profile.  These  integrations  give 


and 


S " 2 ““o 


r.  57T 

C = - — au 
m 32  o 


{L(k)  + I l^dc)}  = I au^P(k)e^“‘' 
{L(k)  + I^Ck)}  au^Q(k)e^‘^'' 


(32) 


(33) 


The  bracketed  terras  represent  the  ratio  of  nonsteady  to  quasi-steady  response 
and  they  will  be  denoted  by  P(k)  for  the  lift  ratio  and  by  Q(k)  for  the 
moment  ratio.  The  integrals  l£  and  I3  are  defined  by 


and 


r 1 _ 1 -jkri 


’ Vn 


1 e dn 


I (k)  = / 1 e ndn 

^ ' Wn 


(34) 


(35) 


Although  these  Integrals  show  some  similarities  to  others  of  hypergeometric 
form,  they  are  not  reducible  to  standard  representations.  They  must  be 
evaluated  numerically. 

The  integrations  for  I2  and  I3  have  been  carried  out  numerically  as 
indicated  in  tlie  Appendix.  The  results  of  this  work  are  shown  in  Figure  5 
which  is  a polar  plot  of  P(k)  and  Q(k)  with  selected  values  of  k showm  on 
the  curves.  Figure  6 compares  P(k)  witli  functions  of  Soars  [7]  and  Horlock 
[5],  referred  to  the  profile  nose  and  with  the  reduced  frequency  based  on 
the  whole  chord.  The  transverse  gust  function  for  cavity  flow  L(k)  [11,12] 
is  also  plotted  in  these  graphs.  In  all  of  these  graphs,  it  is  seen  tliat 
the  steady  state  limit  of  P and  Q as  k >■  0 is  2.  As  pointed  out  previously 
by  Horlock  [5],  this  doubling  is  caused  by  the  fact  that  the  gust,  being  in 
the  direction  of  flight,  affects  the  dynamic  pressure  experienced  by  the 
foil.  But  the  dynamic,  pressure  depends  upon  th.e  square  of  the  velocity  and 
this  fact  is  resnonsible  for  the  limit  noted. 

The  indiclal  admittance  functions  corresponding  to  an  encounter  of  the 
foil  with  a sliarp-cdged  gust  or  unit  step  in  ctiordwisc  disturbance  velocity 
can  be  found  by  direct  superposition  [13, .10],  For  the  lift  response,  one 
can  wri to 


00 
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A^(X)  = f / Re  {P(k)}  dk 


(36) 


For  the  moment  response  one  has 


- f / f5(k)}  dk 


(37) 


Again,  numerical  Integrations  are  required.  An  outline  of  this  work  is  also 
given  in  the  Appendix.  The  results  are  shown  in  Figure  7 which  shov^?s  the 
two  admittance  functions  (36)  and  (37)  plotted  against  X.  The  quantity  X is 
the  distance  traveled  in  chord  lengths  after  the  profile  nose  hits  the  gust. 
Wlien  these  functions  are  known,  the  lift  and  moment  coefficients  are 


and 


C,(X) 


2 “'^o 


i «^>o 


(38) 


(39) 


NONSTEADY  CAVITY-GUST  INTERACTIONS 

It  remains  to  determine  the  effect  of  the  term  w|  = u^e^^^^^  >^/j) 
nonsteady  hydrodynamics  of  the  cavity  flow.  E'rom  the  previous  discussion  of 
the  second-order  boundary  conditions  we  have  seen  that  the  nonstcady  effect 
of  w^  is  ascribable  to  that  part  of  the  complex  acceleration  potential, 
p2(z;2w). 


that  the  boundary  conditions 

for  this  fur 

/ 

2 

w 

0 

0+ 

(i)  Re  {F2}  = - ~ ; x > < 

1 , y = , 

} 

s 

0- 

> / 

< 

(ii)  F2  continuous  at  x = s,  y = 0- 


(iii)  a = a - 0,  at  x = 
X y ’ 


(iv)  V2  = 0 ; 0 < X < s,  y - 0- 
dF 

(v)  Im  {- — } = 0,  0<x<s,  y=0- 
dz  “ - 


We  have  also  observed  tlint  conditions  (iv)  and  (v)  are  somewhat  redundant,  but 
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if  one  knows  that  (v)  is  satisfied  by  F2  he  must  still  be  certain  that  (iv) 
is  explicitly  satisfied.  Finally  it  is  to  be  noted  that  now 


?2 


In  the  following  we  will  Introduce  the  reduced  frequency  k of  Equation  (27) 
and  factor  the  term  w^  into  spatial  and  temporal  parts  so  that 


2 2 „.,x 

''l  “o  "2j^s  2jo)t 

2 2®  ® 


The  temporal  part  can  then  be  elliminated  from  further  consideration  by  writing 
F2  = Fe^j*^*-  and  so  for£h  for  the  other  quantities  defined  above.  The  boundary 
values  which  apply  to  F = $ + iip  in  the  v-plane  are  illustrated  in  Figure  8. 

In  the  2-plane,  the  normal  acceleration  condition  on  the  wetted  surface 
of  the  foil  is  d(p/dx  = 0 and  this  can  be  integrated  along  the  wetted  surface 
to  give 


55  = B(k)  on  y = 0—,  0 < x < s 


Therefore  the  values  illustrated  in  Figure  8 in  the  V-plane  on  the  real 
axis  can  be  replaced  by 


^ O 

u 

^ ^ on  the  upper  and  lower  surfaces  of  cavity 


ip  = B(k) 


on  the  wetted  surface  of  the  profile. 


The  constant  of  integration  B(k)  will  be  determined  later  by  applying  the 
condition  v = 0 on  the  wetted  surface.  For  the  present,  one  observes  that  tlic 
boundary  values  for  F are  now  specified  along  alternate  intervals  of  the 
^-axis  in  terms  of  ? and  iji.  One  can  transform  this  specification  into  a 
boundary  value  problem  on  just  the  real  part  of  an  analytic  function  G(v)  by 
introducing  the  transformation, 


1 ^ ^ _i_  • 

G y -y-  = 1 = + ly 
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For  then,  we  have  on  the  ^ axis: 


c < -1  , g\^yT^  = ? + iij;  ReG  = ; 


for  - 1 < ? < 0 , IG  = $ + lip  ^ ReG  = ’ 


ns: 


for  C > 0 , = $ + lip  ReG  = 


Thus  ReG  is  specified  on  the  entire  real  axis  and  IniG  can  be  found  on  the  real 
axis  by  use  of  the  Hilbert  Transform  [15], 


Substituting  appropriate  values  of  ReG  into  Equation  (A2)  one  finds  that 


-ip/  pf ' /— 

2tt  ^ / V It  I - i IT  J/i-iTTt-c 


O D ; t _dt^ 


' t at 

V t + 1 1 - 


c • 


If  one  substitutes  the  appropriate  signed  value  of  |t|  in  the  above  integrands 
and  also  rearranges  the  limits  of  integration  he  finds,  after  dropping  tlie 
absolute  value  symbol  from  the  dummy  of  integration,  that 


1.C  . - > 1. 

2tt  / 


/ -2jkti  ixr  pK  ‘ 

J t - 1 t + 5 TT  J \/  1 - t t 


+ i p r _‘iL_  . 

2 i , ■ t + 1 t - S 


V.lien  ^ is  on  the  welted  surface  of  the  fell  - 1 < C 5 C ' “ H 
In  this  last  equation  and  find  the  principal  value  of  the  integral  involving 
B(k)  to  be 


1 
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t dt 

1 - t t - n 


The  remaining  integrals  are  no  longer  improper  for  points  on  the  foil  and  so 


ImG(n)  = BOO  + ^ j 


^-2jkt  ^ ^ 

V t-1  t-n 


From  Equation  (41)  one  finds  for  points  on  the  wetted  surface,  with  0 < q < 1, 
that 


$ + i!j;  = 


- ImG  (n)  j - + i j ReG  (n)  I ^ 


Therefore,  the  acceleration  potential  on  the  foil  is 


?(n,0)  - B(k)  - /e 

V h ^ " j J 


'^o  j ;^-2jkt^..,  [~t~  dt  -2jkt^_  r~t  dt 

” \J  ^ ''  ^ ' ~ J V'  t-1  t-n 


. (44) 


» VA  • 

But  the  Kutta  condition  requires  that  ^(l,0;k)  = - e in  accordance  with 
conditions  (1)  and  (ii)  above.  The  addition  or  subtraction  of  a constant,  such 
as  this  one,  to  ^ has  no  effect  on  boundary  condition  (iii)  so  wc  can  write 
Equation  (44)  as 


?(n,o)  = - 


i [I  ‘-'-""S 


I t+1  t+n  / 


r t dt 

V t-1  t-n  f’ 
J i 


(45)  I 


Equation  (45)  can  bo  used  to  calculate  the.  lift  and  moment  on  the  profile.  These 
follow  from  Equation  (17b)  vdiich  is  now  v.’ritLcu  as 


C = - 2 \ 
P 


2-n  / \ t il  t+n 


•2jkt^-^^  ' t i^t_ 

V t-1  t-n 


. (4G) 


In  terms  of  the  present  notation  n = x/s;  and  one  can  write 


c 2 c (n,0)ndn 

k I P 

J 


and 


C = 
m 
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(47) 


In  the  second  of  equations  (47)  the  moment  is  measured  about  the  nose  of  the 
profile  and  is  taken  as  positive  in  the  direction  of  a increasing.  The  use 
of  Equation  (46)  in  Equations  (47)  leads  to  a number  of  integrations  which 
can  be  obtained  in  closed  form.  For  example,  those  involving  B(k)  are  among 
these.  In  other  instances  one  interchanges  the  order  of  integration  in 
Equations  (47)  and  finally  obtains 


o M -2jk' 

^ n , 2 1 1-e  1,0  / 

Cj»  jB(k)  -yt-j 


and 


i 


-2jkt^ 


^,/t(t+l)  - 


r 2f  i 


dt  > 


C = - ~ b(k)  - u^ 
m 32  o 


,-2jk 


4k 


2jk 


16  / t+1 


/t^  + (f  “ t“)  ; t(t+l)  ; dt 


I-'  f-L- 

1 / t+1 


dt 


1 \ 


(48) 


(49) 


It  now  remains  to  evaluate  the  coefficient  B(k)  by  applying  condition  (iv) 
above.  As  a preliminary  to  tliis  evaluation  we  note  tliat  liaving  obtained 
Equation  (45)  we  now  know  ?i(n,0)  cve.rj’where  on  the  real  V-axis.  Moreover, 
ip(n,0)  = B(k)  on  the  wotted  surface,  0 < f|  < 1,  but  we  have  not  as  yet  found 
it  on  other  parts  of^the  real  V-axis,  However  for  our  purposes  it  is 
necessary  to  find  ImF(V)  on  the  upper  half  of  the  Imaginary  v-axis.  This 
deterniinatlon  is  facilitated  if  F(v)  is  treated  as  the  sum  of  two  functions. 


F(V)  = F(v)  + H(v) 


where 


ReF(f|,0)  = 0,^<-l,5>0  (on  tlie  cavity) 

ImF(^,0)  ==  B(i:)  > “3  f f 0 (on  the.  foil) 


J 
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and  where 


ReH(5,0)  = - e , ^ < -1  , ^ > 0 (on  the  cavity)  , 

ImH(£;,0)  = 0 , -1  < C < 0 (on  the  foil) 


(50) 


The  boundary  value  problem  for  F has  the  same  form  as  the  first  order  steady 
flow  problem  treated  previously.  Therefore  we  can  write  dow'n  the  analytic 
function  F(?)  without  further  ado.  Thus  in  the  ^ and  the  V-planes  we  have 


F = 


-2iB 

S-1 


iB  1 - \ / 1 + 


Then  if  one  puts  v = iri  for  points  on  the  upper  half  of  the  imaginary  axis  of 
the  V-plane,  he  can  write 


if'  = - B(k) 


IniH(0,n) 


(51) 


In  order  to  find  the  function  II  wc  employ  the  transformation  (41)  and 
the  Hilbert  transform  (42)  to  determine  Imll  on  the  real  axis  off  the  foal  in 
the  V-planc.  Then  Imll  is  known  everywhere  on  the  real  axis  and  its  value  on 
the  upper  half  of  the  imaginary  axis  in  the  v-plane  can  be  obtained  with  the  aid 
of  the  Poisson  integral.  In  particular  we  find  from  (41)  and  (50)  that 


' "L/  -2jkrd  _ , 
T V ;; — ' ■ 5 < -1  ; 


Id 


ReC  =--< 


-1  < C < 0 ; 


' ^ r T -1  r2 

; - , 0 < 5 . 

I ^ •'  C + 1 


Then  dt  follows  ^or  t,  =•  - q < -1  that 

i OO  00 

^ ' ri  ^ r~r  -2jkt^  f 

'./a  - dt  i ' ■ 
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-2Jkt 


2ti  ^ q I ( V t-.1  q-t 
1 ‘ 

I 


J 


7 t+3  q-lt 


dt 


Note  that  when  f,  = -1  (q'l)  IniH  = 0.  Thus  Ihe  Kul  I a condition  is  satisfied. 
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= i-,  /5T 1 at . P Aa:  dt 

2tt  VC  I / V'  t-1  C + t J V t+1  t - c 


Now  that  ImH  is  specified  all  along  the  real  axis  we  can  use  the  Poisson 
integral, 


ImlKO.ln)  - 5 dt  , 


to  write 


11.11(0.1,1)  " A '■'1  ’’  f 'V'H^  I 

2^  J ^ J 


J V q (q-t)(q^+n  ) 


_ ,^t  -2ikt^  ' ••  q-1  <^qdt 

' / e i \ / 


2 2 I t-1 

(tJ-q)(qAn  ) J ’ 


„-2jkt  I / C+1  dCdt  

J ^ ^ (t+c)(C^+n^ 


T-,  'it_  -2jkt^  r /|k  dgdt 

j j V (t-c)(f/+n^) 


Therefore  it  folJows  from  Equation  (51)  that  the  desired  result  is 


iP(0,in)  = - B(k)  1 V 1 ij  + ImH(0,in)  . 


One  can  transform  this  expression  from  the  upper  half  of  the  imaginary  v-axis^ 
to  the  negative  ) cal  £o;is  in  tlie.  z-p]ano  by  means  of  the  relation  (x/s)  - “ q . 
The  arguments  leading  to  Equation  (29)  apply  here  too  so  that  condition  (iv) 
of  zero  downwash  V2  on  the  foil  becomes 


0 = - B(k)  + 2jk  B(k)  f i 1 - ; " J Vl  + -i  + 1 j dt  - 


- 7 -i  V ^ 

2jl:  ! c H(t,0)clt 
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As  before  we  have  changed  the  variable,  of  integration  from  - (x/s)  to  t and 
the  limits  of  integration  have  a] so  been  altered  to  suit  this  change.  The 
function  H(t,0)  is  as  defined  in  Equation  (52)  with  it)  (-x/s)  ->■  t also. 
Evidentially  the  function  B(k)  is  defined  by 


-2jk  J 


e"^^‘"‘^H(t,0)dt 


B(k) 


1 + j(2k)I^(2k) 


(53) 


where  the  Integral  l2(2k)  lias  been  defined  previously  by  Equation  (30)  or  by 
L(2k)  in  Equation  (31).  As  it  stands  the  integration  in  the  numerator  of 
Equation  (53)  involves  triple  Integrals.  However  if  one  inverts  the  order  of 
integration,  the  last  integration  can  be  written  as 


CO 


j 


7 2 

[-Ci(2k^“)  + j Si(2kC  ) 


Therefore  Equation  (53)  can  be  written  as 


B(k)  = -2jk  ~jL(2k) 
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Equation  (54)  can  be  used  1ti  Equat.i.ons  (^‘8)  and  (49)  to  complete  tlie  determination 
of  the  non.'Ueady  lift  and  moment  resulting  from  tiie  gust-cavity  interaction 
for  sinusoidal  gu.sts.  Once  t'.iesc  quantities  have  been  evaluated  numerically  tlie 
corr  es.pond  tnr,  rcspon.ses  to  a shaip-ed{;ed  streawwise  gust:  may  also  be  d<!tcrmincd 
for  Ibis  jmrt  of  ihe  solntir>n  as  v.'as  done  before  for  the  other  i>art  . The  se 
solutions  ean  then  be  coiiiln'ned  vi'Ctorially  to  give  tlie  total  nonsteady  response 
in  term.s  of  I and  (he  ratio  a/u(,  for  sinusoidal  girOs.  I'or  siiarp-edged  gusts, 
curves  similar  to  I'if.ure  7 ean  he  obtained  for  various  values  of 
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CONClAlSIOr^S  AND  RFCOMMEN'DATIONS 

The  present  study  has  shown  that  there  is  a cavity  flow  analog  of  the 
Horlock  function  of  airfoil  theory  which  gives  the  response  of  the  profile 
to  strearawise  sinusoidal  gusts.  This  nonsteady  response  function  for  flov;s 
at  zero  cavitation  number  has  been  evaluated  numerically  for  the  entire  range 
of  reduced  frequencies  for  both  lift  and  moment  on  the  profile  and  these 
results  have  been  used  to  derive  the  transient  lift  and  moment  on  the  profile 
after  it  encounters  a sharp  edged  gust.  It  is  found  that  the  results  computed 
can  be  ascribed  mainly  to  the  direct  action  of  the  streamwise  velocity 
fluctuations  upon  the.  inclined  wetted  surface  of  the  hydrofoil.  This  is  tiie 
only  effect  found  in  airfoil  theory. 

It  has  been  found  theoretically  that  the  cavity  flow  causes  the  hydrofoil 
to  experience  an  added  nonsteady  effect  which  is  not  found  in  airfoil  theory. 
This  added  reaction  occurs  at  twice  the.  frequency  of  the  input  frequency  of 
the  gust.  This  part  of  the  solution  has  been  determined  explicitly  in  the 
present  study.  It  is  independent  of  attack  angle  whereas  the  part  of  the 
solution  which  was  determined  first  depends  on  both  gust  amplitude  and  attack 
angle.  Thus  one  can  assign  tlie  cause  of  this  nev/  nonsLeady  effect  to  be 
primarily  an  interaction  between  the  cavity  and  the  gust.  To  date  numerical 
analysis  of  this  part  of  tlie  tlieory  has  not  been  carried  out. 

Tlie  formulae  v.’hlch  have  been  derived  for  this  second  part  of  the  solution 
indicate  that  the  numerical  work  for  its  evaluation  will  be  rather  elaborate. 
Once  this  work  has  been  done  the  total  response  of  the  profile  can  be  obtained 
as  a vectorial  combination  of  the  two  sets  of  data  for  various  ratios  of  steady 
attack  angle,  to  gust  amplitude  over  the  range  of  reduced  frequencies.  For  the 
total  response  to  the  sharp-edged  gust,  results  for  the  two  parts  may  be  added. 

The  ability  of  an  oscillating  surface  to  extract  energy  from  the  flow 
depends  critically  upon  the  pliasc  relationships  between  tlie  nonsteady  reactions 
and  the  oscillatory  input,  T'ne  amplitude  of  these  nonstcady  reactions  need 
not  be  large  because  it  is  the  stability  of  the  response  which  is  important. 
Therefore,  it  is  essential  to  evaluate  this  second  part  of  the  sol ution  in 
order  to  explore  this  important  aspect  of  this  new  theory.  It  is  conceivable 
that  the  effects  discovered  during  the  course  of  this  study  could  be  important 
for  analyses  of  the  structural  dynamics  of  fully  cavitating  hydrofoils  in  a 
seaway.  For  these  reasons  it  is  recommended  that  upon  completion  of  the  present 
study  the  theory  be  extended  to  include  free  surface  effects. 


APPENDIX 

NIIMKRICAL  CALCULATION, S 

As  long  as  t he  reduced  frequency  is  not  too  large,  numerical  eval  u.at.ions 
for  I2  and  I3  as  defined  by  liquations  (.lA)  and  (35),  can  proceed  by  standard 
methods.  For  the  present  Simpscui's  rule  int  egratii'iis,  l;-valiics  ranged  frinn 
zero  to  twenty-five.  This  rather  large  range  \;js  scloctcnl  in  order  to  permit 
the.  most  efrcctivu  joining  v.’itli  results  obtained  by  means  of  asymptotic 
expansions  for  riduced  frequencies  rang.ing  from  ten  to  infinity.  It  also 
permits  a check  (’11  the  validity  of  tlie  numerical  integrn Li ous . 


It  was  found  tiiat  ror.pulcr  Lime  could  tie  jcduccJ  if  oru  uses  t lie  fact  llint 
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~ “ 5ir/64,  for  then  one  can  write 


I,  -?+  / -\/4=-  i 


Re  I^ 

= 7 + 

' j 

0 
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= 5tt  ^ 

Re  I^ 

64 

(coskri  - 1)  dr)  , Im  I, 


V n 


/ VVn 


slnkr)dri  , 


I3  = U + ; - - fii)  (cosk  - 1)  dn  , Im  I3  = - / T\'^-  ~ Yn)sinkndn  . 


3/2 


j 


Moreover,  because  all  four  integrands  above  have  branch  points  at  both  r)  = 0 
and  n = 1»  ii  was  found  lielpful  to  devise  special  finite  difference  integration 
formulas  for  use  near  these  points.  All  integrations  were  carried  out  witli 
extended  precision.  They  also  used  an  accuracy  test  to  determine  an  appropriate. 
Interval  size,  for  each  k-value  and  for  each  of  the  fo\ir  integ.rals.  In  order  to 
obtain  values  of  P and  Q,  one  must  also  have  values  of  the  function  h(k).  'J’hese 
numerical,  data  v;ere  taken  from  the  tabulation  of  Reference  12  and  used  in  the 
expressions 


and 


p = hOO  + 4 i,(k) 
ir  2 

Q - Uk)  + f 13(10  . 


For  values  of  k ranging  from  ten  to  infinity,  the  method  of  steepest 
desce.nt  [16,17]  was  used  to  find  asymptotic  expansions  for  I2  and  1 3.  As  one 
might  ex[>ect,  a contribution  to  the  value  of  these  lnt''„rals  is  obtained  from 
both  branch  points  of  the  integrand.  The  path  of  integration  ns  specified  in 
terms  of  the  real  variable,  ri  , between  0 and  1,  corresponds  to  a level  curve. 
The  results  obtained  by  this  procedure  are: 
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r(7/2) 

(jk) 


In  order  to  calculate  the  response  function  P(k)  and  Q(k),  the  preceding 
results  must  be  supplemented  by  asymptotic  values  of  the  function  L(k)  as 
indicated  in  the  above  Kquations  for  P and  Q.  3’hese  additional  values  can 
be  calculated  directly  from  tlie  asymptotic  expansion  of  the  quantity 
1 + jklj^  [13].  In  order  to  obtain  values  which  provided  adequate  accuracy, 
it  seemed  necessary  to  obtain  more  ternis  from  the  expansion  of  1 + jkli  than 
is  provided  in  Reference  13.  The.  required  extension  of  the  formula  is 


r(5/4)  1 r(7/4) 

,1/4  8 ,3/4 

(jk)  (jk) 


1 + jkl^(k)  ~ -^<r(3/4)  (jk)^''^  + I 


r(9/4) 
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(jk) 


5 rqi/4) 

128  ,7/4" 
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An  additional  term,  , was  also  found  in  order  to  determine  the 

best  number  of  ternis  to  retain  in  the.  expansion.  It  was  found  that  this 
additional  term  x,rould  have  reduced  the.  accuracy  of  the  calculations  and  that, 
for  values  of  reduced  frequency  exceeding  ten,  little,  if  any,  difference 
is  found  betv.'een  L(k)  values  calculated  from  the  expansion  of  1 H-  jklj^,  as 
shown  above  or  with  only  the.  first  four  terms  as  originally  given  in 
Reference  13.  Hov.'ever,  it  was  found  helpful  to  recompute  all  the  coeffi- 
cients of  k in  this  formula  to  more  decimal  places, 

A tabulation  of  all  numerical  results  for  the  sinusoidal  gust  response 
is  shown  in  Table.  1.  These  values  of  the  complex  functions  defined  by 
Equatlotis  for  P(k)  and  Q(k)  are  sliown  plotted  against  reduced  frequency  in 
Figures  8 and  9.  It  should  be  recalled  in  connection  wi tli  these  results  that 
all  moments  and  forces  arc  referred  to  the  nose  of  tlie  foil  and  that  the 
reduced  frequency  is  based  on  the  profile  chord  and  not  upon  the  half  chord 
as  is  common  in  airfoil  theory. 

TJie  indlcial  admittance  functions  corresponding  to  an  encounter  of  tlie 
foil  with  a sharp-edged  gust  or  unit  stop  in  chordwise  disturbance  velocity 
can  be  found  by  direct  superposition  [13,14],  For  the  lift  response,  one  can 
write 
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A (A) 
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dk 


(39) 
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In  these  formulae  the  quantity  X is  the  distance  traveled  in  cliords  after  tlie 
profile  nose  hits  the  gust.  We  have  already  noted  that  Re  {P}  ->■  2 and 
Re  {Q}  ->  2 as  k -»  0.  Moreover,  the  behavior  of  both  Re  {P}  and  Rc  {q}  is 
dominated  by  Re  {h(k)}  as  k ->■  It  has  been  found  that 
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Evidentally,  Re  {P}  and  Re  {Q}  behave  well  enougli  at  infinity  so  that  the 

Fourier  transforms  for  A„  and  A converge  and  one  can  v/rite  them  as: 
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2 - Re  Q 


This  transformation  takes  care  of  the  integrands  at  the  origin  so  tliat  asymptoti c 
representations  for  these  integrals  as  X are  easily  obtained.  It  is 
necessary  to  approximate  tlie  quantities  2 - Re  P and  2 - Re  Q near  k - 0 by 
suitabie  polynomials.  From  Table  1,  it  is  found  that 


2 - Rc  P = .20501; 


2 - Rc  Q = .2025k  + .12501; 


Thcfje  forms  appear  to  be  consistent,  with  the  accuracy  of  tlie  tabulated  values. 
Now,  it  is  knovm  [18]  tliat  if  F(x)  and  all  of  its  derivatives  exist  as  ordinary 
functions  for  x ^ 0,  and  are  well  behaved  at  infinity,  then  the  asymptotic 
cxqTansion  of  the  Fourier  transform  of  F is 
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FfOl  F"(0)  , F (0) 

F(x)  sin  2Tf  xy  dx  ~ n + - --V-  - ... 

^ ^ (2ty)^  (2iiy)-^ 


In  order  to  apply  this  result  to  the  present  problem,  we  put 


F(0)  -=  ^ .2050  or  F(0)  = - .2025 
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and  X = 2ny.  Thus,  one  has 


These  expansions  are  used  to  suppiement  numerical  integrations  i.’hich  are  based 
upon  the  values  of  Table  I and  the  asymptotic  expansions 


Results  of  such  calculations  for  Aj^CX)  and  Ajj^(X)  arc  presented  In  Table  II. 
Graphs  of  Aj^  and  are.  shown  in  Figure  7. 
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FIGURE  1 - Flow  Geometry  and  Physical  Quantities 


• I 

FIGURE  2 - Mathematical  Representation  of  the  Slender  Body 
Approximation  of  the  Cavity  F.lov.' 


- 37  - 


Juno  4,  1974 

BRPrlhm  ^ 

j 


I 


FIGURE  3 - Foil  and  Cavity  in  the  V--I'lana 
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FIGURE  4 - Profile  and  Cavity  in  tlic  C''Flanc 
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A second-oi'der  linearised  theory  of  two-diTT'ensi.onal  cavity  flows  is 
used  to  study  the  nonsteady  response  of  hydrofoil  lift  and  morent  tc'  sinusoidal 
and  sharp-edp,ed  streamvi.se  pusts.  -ylli e?~c-  p.ust  patterns  are  assumed  to  be  frozen 
and  to  be  convecLed  at  the  steady  frec'-stronin  veloci t?7~'3 Tlie  iiiialysis  is  restricted 
to  the  case  of  Z6!ro  cavitation  number  cor  respond  i n;;  to  an  infinitely  long  cavity 
in  an  oU1cn^’ise  unbo^'iidcd  flow.  Po}  sinusoidal  gusts,  lift  and  moment  response 
functions  are  presented  for  tb.e  c.niiTo  range  of  reduced  frequency.  In  addition, 
transient  lift  and  momeiU  responses  aic  t.cbiilatrci  for  tlie  reactions  whidi  occur 
after  the  foil  encounters  a sliar]i-eJg,i_d  g.ust.  These  calculations  are  carried 
out  for  those  terms  in  the  solution  vlri ch  result  from  the  nonsteady  dovnwasli  on 
tlic  wetted  surface  of  tb.e  foil  which  are  due  to  t!ie  dii'cct  action  of  (be  gust 
on  the  inclined  vetted  surf.ace.  They  provide  a direct  cavity  flow  analog  of 
the  Horlock  function  of  airfoil  theory,  'ilie  present  study  of  a cavity  flow  has 
revealed  an  added  i.ou.steady  effect  r«.T.u1tirg  f roir.  an  interaction  between  the 
gust  and  the  cavii.y.  Thcorc tica.l  rc.sults  for  thi.c  part  of  the  solution  are 
presented  in  this  re^rt  win" ch  show  a response  at  twice  t're  input  frequency 
of  a sinusoidal  g .s t . \Nur.erical  evaluation  of  this  part  of  tr;e  theory  has  not 
Leon  carried  ou(. . 
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